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Abstract 

We report a first principles systematic study of atomic, electronic, and magnetic properties 
of hydrogen saturated silicon nanowires (H-SiNW) which are doped by transition metal (TM) 
atoms placed at various interstitial sites. Our results obtained within the conventional GGA+U 
approach have been confirmed using an hybrid functional. In order to reveal the surface effects 
we examined three different possible facets of H-SiNW along [001] direction with a diameter 
of '^2nm. The energetics of doping and resulting electronic and magnetic properties are ex- 
amined for all alternative configurations. We found that except Ti, the resulting systems have 
magnetic ground state with a varying magnetic moment. While H-SiNWs are initially non- 
magnetic semiconductor, they generally become ferromagnetic metal upon TM doping. Even 
they posses half-metallic behavior for specific cases. Our results suggest that H-SiNWs can be 
functionalized by TM impurities which would lead to new electronic and spintronic devices at 
nanoscale. 

Introduction 

Among the research for building blocks in nano-devices, silicon nanowires (SiNW) are attracting 
increasing interest due to their remarkable physical, electronic, thermal, and chemical proper- 
ties.'^ Currently, rodlike, oxidation resistant SiNWs can be fabricated with a diameter down to 
1 nm where quantum confinement effects are also observable.'^ The compatibility with current 
silicon based technology*^ makes SiNWs even more attractive and they can enter various device 
applications such as field effect transistors,'^ light emitting diodes,^ lasers,'^ nanosensors,'^ etc. 
The advances in synthesis will possibly increase the number of potential applications into various 
other fields. 

At the theoretical level, bare and hydrogen saturated silicon nanowires (H-SiNW) have been ex- 
tensively analyzed by using first principles calculations .^^^^ It is shown that while bare SiNWs are 
in general metallic, they become semiconductor (insulator) when saturated with hydrogen.'^^EHIl 
The electronic band gap can be engineered by varying diameter, growth direction and cross sec- 
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tion due to the quantum confinement effects.'^'^^' More recently, superlattice structures leading to 
confined states induced by diameter modulated SiNWs'^and merged Si and Ge nanowiresl^have 
been predicted. 

Furthermore, growing research interest has also been devoted to functionalization of SiNWs 
with various dopants in order to study the chemical and biological sensitivities J 1 1 1 ^ 1 1 22 ] j^^g ^j^^ 
p-doped H-SiNW can provide excess carriers required in device applications like diodes and tran- 
sistors . ^^^^1 Recently, Wu et al.'^have reported the room temperature magnetism of Mn-implanted 
SiNWs which opens the field to spintronic applications through transition metal (TM) doping. The- 
oretically Mn impurities in [111] SiNWs and also external adsorption of various TM atoms on 
H-SiNWs are reported.!™! 

In this paper we present an extensive first principle study on hydrogen saturated silicon nanowires 
which are doped by transition metal atoms (Ti, V, Cr, Mn, Fe, and Co) placed at various interstitial 
sites. The energetics of doping and resulting electronic and magnetic properties are examined for 
all alternative configurations. In order to reveal surface effects, three different facets along [001] 
direction are considered. The magnetic ground state is determined by considering non-magnetic, 
ferromagnetic, and anti-ferromagnetic configurations. We also considered calculations by using a 
new type of hybrid functional which in turn are compared with the results obtained by GGA-i-U. 

Methodology 

We have performed first-principles plane wave calculations within density functional theory(DFT)'^ 
implemented in VASP code.*^ All calculations for non-magnetic, ferromagnetic and anti-ferromagnetic 
states are carried out by using the projector augmented wave (PAW) potentials.'^^'^The exchange 
correlation potential has been approximated by generalized gradient approximation (GGA)^ 

To model our silicon nanowires, we have used a supercell approach, including a 10 A vacuum 
space between wire replica along the two directions perpendicular to the wire axis. For doped 
nanowires, the impurity TM atom is periodically repeated along the [001] direction corresponding 
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to the wire axis. In most calculations, one TM atom is introduced per basic wire unit, corre- 
sponding to a distance between consecutive impurities of 5.4 A. Even we are not aiming to obtain 
impurity levels, in order to test the size dependence of our results, for specific cases we doubled 
our supercell size along [001] direction, so doubling the distance between an impurity TM atoms. 
The results were reasonably similar,!^ confirming the main trends and our conclusions. As ex- 
pected, Eb slightly increases with lengthening the supercell size owing to the reduction in TM-TM 
interaction (which enhances TM-SiNW interaction), but electronic band structure profiles are not 
affected. 

In the self-consistent potential and total energy calculations the Brillouin zone of supercell is 
sampled in the k-space within Monkhorst-Pack scheme'^by (1x1x15) mesh points. A plane-wave 
basis set with a kinetic energy of 500 eV has been used. All atomic positions and lattice constant 
along the wire axis are optimized by using the conjugate gradient method where total energy and 
atomic forces are minimized. The convergence in energy is chosen as 10^^ eV between two ionic 
steps, and the maximum force allowed on each atom is 0.01 eV/A. 

A second set of calculations has been performed using the so-called Bl-WC hybrid func- 
tional,^ which mixes 16% of Hartree-Fock exchange with Wu-Cohen GGA^^ within Bl scheme 
These calculations have been performed using the linear combination of atomic orbitals method, 
as implemented in the CRYSTAL code.^^ The Bl-WC hybrid functional used in this work was de- 
veloped earlierl^ and it was tested for H-SiNWs.*^ We used localized Gaussian-type basis sets 
including polarization orbitals and considered all the electrons for Si^ H, and Fe.l^ The 

Brillouin zone integrations were performed using a 1x1x15 mesh of k points and a secondary 
2x2x30 Gilat k-mesh was used in the calculation of Fermi energy and density matrix. The self- 
consistent-field calculations were converged until the energy changes between interactions were 
smaller than 10^^ Hartree using a Fermi smearing of 0.00032 Hartree (~100 K). An extra-large 
predefined pruned grid consisting of 75 radial points and 974 angular points was used for the 
numerical integration of charge density. The level of accuracy in evaluating the Coulomb and ex- 
change series is controlled by five parameters: the values used in our calculations are 7, 7, 7, 7, 
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and 14. The Bl-WC calculations have been performed for the optimized structures obtained with 
the VASP code. 



H-SiNW{45) H-SiNW(57) H-SiNW(61) 




Figure 1 : (Color online) The top and side view of hydrogen saturated silicon nanowires along 
[001] direction with different facets for N=45,57 and 61 where N corresponds to number of silicon 
atoms in the unitcell of nanowire. 



Hydrogen Saturated Silicon Nanowires 

Before considering interstitial TM doping, we studied structural and electronic properties of pro- 
totype bare [SiNW(N)] and hydrogen saturated [H-SiNW(N)] silicon nanowires for N=45, 57, 61 
where N is the number of silicon atoms in the unit cell of nanowires ([T]). These H-SiNWs corre- 
spond to three possible facets along [001] direction and their diameter changes between 1.5-2.1 
nm.'^ As already discussed in the introduction part, H-SiNWs as small as 1 nm have been fabri- 
cated^ and at such small sizes edge effects become more important.l^^^^ 

The initial atomic positions of SiNWs are taken from silicon bulk crystal. Depending on the 
cross-section three alternative facets are possible along [001] direction. Upon ionic relaxation 
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H-SiNW{45) H-SiNW(57) H-SiNW{61) 




Figure 2: (Color online) The band structures of hydrogen saturated silicon nanowires along [001] 
direction for N=45,57, and 61. The electronic energy band structures and the gaps (Eq) are calcu- 
lated using GGA. The values given in parenthesis are obtained using Bl-WC calculations. 

by minimizing both the total energy and forces on the atoms, the ground state configurations are 
obtained. The cohesive energy (Ec) of the SiNWs are calculated by the definition Ec=ET[Si]- 
ET[SiNW(N)/N] where ExLSi] and ET[SiNW(N)] is the total energy of free silicon atom and 
SiNW(N), respectively. Ec is calculated as 4.91, 4.97, and 4.99 eV for N=45, 57, and 61, respec- 
tively. Interestingly, SiNW(61) which has sharp corners has the highest Ec which is in agreement 
with the results of Cao and coworkers. Clearly, the presence of edges has a profound effect on the 
surface reconstruction of SiNW and thereby its electronic structure and stability .'^^'^ Accordingly, 
while thick SiNWs prefers cylindrical or prism shape with a core that preserves diamond structure, 
the cross section with sharp comers becomes more favorable for very thin nanowires.!^ 

As a next step, we saturate the dangling bonds of silicon atoms on the surface with hydro- 
ggj^ 151171 j(- |g observed that while SiNWs are all metallic due to the surface states, they become 
semiconductors upon hydrogen saturation as shown in |2] The bang gap (Eg) tends to decrease 
from N=45 (Eg=1.79 eV) to N=61 (Eg=1.16 eV) due to an increase in diameter (despite the cross 
sections are different). This is consistent with predicted quantum confinement effects. 
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It is well known that while calculations based on DFT typically yield very accurate structural 
properties for silicon and other semiconductors, the standard approximations (LDA or GGA) lead 
to a significant underestimation of the electronic band gaps due to intrinsic failure in handling self- 
interactions. ^21 Accordingly we recalculated the band gaps by using Bl-WC hybrid functional,'^ 
the details of which are given in the methodology part. We noticed that the band structure pro- 
file that are obtained by GGA and Bl-WC are very similar, and the correction results in a rigid 
shift of conduction and valence bands corresponding to a significant increase in the band gap. 
The corrected Eq's are 2.53, 2.38 and 1.81 eV for N=45,57, and 61, respectively as shown in|2j 
The obtained values are in agreement with the experimentally available datal^ and theoretical GW 
corrected results'^ for H-SiNWs within same diameter range.l^ 

Interstitial Transition Metal Doping 

Previously, we showed that TM atoms can strongly bind to H-SiNW(N) surfaces without deform- 
ing the nanowire structure .'^^'^ As the doping occurs at high temperatures there is a strong possibil- 
ity that TM atom can also diffuse to interior regions. Accordingly we extend our previous analysis 
and consider various possible interstitial sites for six TM atoms (Ti, V, Cr, Mn, Fe, and Co) and an- 
alyze energetics, structural, electronic, and magnetic properties of TM doped H-SiNWs, labeled as 
H-SiNW(N)-i-TM(s), where TM is the type of transition metal atom and s is the interstitial doping 
site as shown in [3} The obtained results, namely bond distances (dsi-xivi). binding energies (Eb), 
magnetic moments(/i) and the differences (AE) between the total energies of spin-unpolarized (su) 
and spin-polarized (sp) states are summarized in[T]and[2j 

Structure and Energetics 

As a first step, we consider different positions from subsurface to core regions at which TM atoms 
can settle. The projection of these sites are shown in|3j This corresponds to one impurity atom for 
6 layers of H-SiNW ([T]). For each site there are two possible geometries for TM doping: the center 
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(a) H-SiNW{45)+TM H-SiNW{57)+TM H-SiNW{61)+TM 




Figure 3: (Color online) (a)The possible interstitial adsorption sites of transition metal (red 
spheres) atoms labeled 1-5 in H-SiNW(N) for N=45,57 and 61 and (b) two alternative doping 
geometries. 



Table 1: The bond distances, binding energy, magnetic moment per cell, total energy differ- 
ence of polarized and unpolarized states for subsurface (s=l and 2) sites. 





s 


N=45 


N=57 


N=61 






d(A) 


Eb(eV) 




AE 


d(A) 


Eb(eV) 




AE(eV) 


d(A) 


Eb(eV) 




AE 


Ti 


1 


2.38-2.85 


3.42 


0.00 


0.00 


2.38-2.86 


3.59 


0.00 


0.00 


2.47-2.98 


4.34 


0.00 


0.00 




2 


2.43-2.84 


3.59 


0.00 


0.00 


2.45-2.80 


3.92 


0.00 


0.00 


2.42-2.85 


4.35 


0.00 


0.00 


V 


1 


2.34-2.77 


2.96 


1.01 


0.02 


2.34-2.78 


3.07 


0.83 


0.02 


2.45-2.83 


3.82 


2.97 


0.02 




2 


2.42-2.78 


3.09 


1.21 


0.02 


2.39-2.75 


3.38 


1.01 


0.01 


2.39-2.78 


3.39 


1.39 


0.02 


Cr 


1 


2.29-2.73 


1.42 


3.75 


0.09 


2.31-2.97 


1.72 


3.17 


0.21 


2.42-2.80 


2.46 


2.94 


0.46 




2 


2.45-2.75 


1.54 


3.72 


0.20 


2.43-2.73 


1.87 


3.51 


0.26 


2.41-2.81 


1.89 


3.46 


0.27 


Mn 


1 


2.27-2.77 


1.75 


2.53 


0.05 


2.27-2.86 


1.89 


2.44 


0.10 


2.40-2.77 


2.66 


2.39 


0.42 




2 


2.42-2.77 


1.93 


3.00 


0.35 


2.41-2.75 


2.21 


3.00 


0.42 


2.39-2.81 


2.22 


3.00 


0.39 


Fe 


1 


2.23-2.74 


3.18 


1.55 


0.19 


2.22-2.85 


3.28 


1.46 


0.19 


2.35-2.74 


3.83 


1.39 


0.35 




2 


2.39-2.73 


3.26 


2.00 


0.35 


2.37-2.73 


3.45 


2.00 


0.31 


2.36-2.81 


3.44 


2.00 


0.28 


Co 


1 


2.17-3.01 


3.94 


0.00 


0.00 


2.19-2.97 


4.16 


0.00 


0.00 


2.30-2.96 


4.76 


0.00 


0.00 




2 


2.39-2.70 


3.80 


1.00 


0.08 


2.37-2.77 


3.99 


1.00 


0.03 


2.37-2.89 


4.02 


1.00 


0.02 
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Table 2: The bond distances, binding energy, magnetic moment per cell, total energy differ- 
ence of polarized and unpolarized states for core (s=4 and 5) sites. Results indicated by 
correspond to hexagonal symmetry (See the text). 





s 


N=45 


N=57 


N=61 






d(A) 


Eb(eV) 




AE 


d(A) 


Eb(eV) 




AE(eV) 


d(A) 


Eb(eV) 




AE 


Ti 


3 


2.45-2.78 


3.65 


0.00 


0.00 


2.44-2.84 


3.86 


0.00 


0.00 


2.42-2.85 


3.95 


0.00 


0.00 




4 


2.45-2.78 


3.65 


0.00 


0.00 


2.47-2.80 


3.75 


0.00 


0.00 


2.45-2.79 


3.72 


0.00 


0.00 




5 


2.45-2.83 


3.72 


0.00 


0.00 


2.48-2.73 


3.57 


0.00 


0.00 


2.47-2.79 


3.78 


0.00 


0.00 


V 


3 


2.44-2.73 


3.20 


2.52 


0.08 


2.42-2.80 


3.33 


1.36 


0.03 


2.39-2.78 


3.39 


1.40 


0.02 




4 


2.44-2.73 


3.20 


2.52 


0.08 


2.43-2.76 


3.24 


1.67 


0.03 


2.42-2.77 


3.28 


1.32 


0.09 




5 


2.41-2.79 


3.19 


1.41 


0.02 


2.46-2.79 


3.19 


2.42 


0.12 


2.40-2.80 


3.21 


2.05 


0.03 


Cr 


3 


2.45-2.75 


1.78 


4.00 


0.40 


2.42-2.79 


1.78 


3.61 


0.25 


2.41-2.81 


1.89 


3.46 


0.27 




4 


2.45-2.75 


1.78 


4.00 


0.40 


2.45-2.76 


1.79 


3.82 


0.34 


*2.44-2.80 


*1.83 


*4.00 


*0.34 




5 


2.43-2.78 


1.73 


4.00 


0.33 


2.47-2.70 


1.77 


4.00 


0.38 


2.32-2.46 


1.91 


4.00 


0.38 


Mn 


3 


*2.37-2.39 


*1.45 


*3.00 


*0.62 


2.40-2.79 


2.18 


3.00 


0.47 


2.39-2.81 


2.22 


3.00 


0.39 




4 


*2.32-2.38 


*1.54 


*3.00 


*0.60 


2.42-2.77 


2.14 


3.00 


0.45 


*2.33-2.37 


*1.49 


*3.00 


*0.35 




5 


2.40-2.79 


2.09 


3.00 


0.52 


2.44-2.69 


2.06 


3.00 


0.36 


*2.31-2.38 


*1.61 


*3.00 


*0.32 


Fe 


3 


*2.26-2.38 


*2.74 


*2.00 


*0.44 


2.38-2.77 


3.41 


2.00 


0.38 


2.36-2.81 


3.44 


2.00 


0.28 




4 


*2.26-2.39 


*2.84 


*2.00 


*0.40 


2.39-2.76 


3.38 


2.00 


0.34 


*2.27-2.37 


*2.81 


*2.00 


*0.41 




5 


2.38-2.78 


3.34 


2.00 


0.45 


2.42-2.69 


3.35 


2.00 


0.27 


*2.26-2.39 


*2.89 


*2.00 


*0.29 


Co 


3 


*2.28-2.34 


*3.60 


*1.00 


*0.06 


2.36-2.79 


3.92 


1.00 


0.07 


2.37-2.89 


4.02 


1.00 


0.03 




4 


*2.28-2.33 


*3.68 


*1.00 


*0.02 


2.37-2.82 


3.94 


1.00 


0.06 


2.37-2.89 


3.64 


1.00 


0.04 




5 


2.37-2.75 


3.85 


1.00 


0.08 


2.42-2.70 


3.95 


1.00 


0.05 


*2.28-2.34 


*3.83 


*1.00 


*0.03 



of a cage where TM atom is surrounded by ten silicon atoms and the center of hexagon where TM 
atom is in the middle of six silicon atoms forming an almost planar geometry as shown in |3] For 
all the cases the cage geometry yields stronger Eb due to an increase in coordination number of 
TM atom. For most of the cases the hexagon geometry is not stable and TM atom moves to center 
of silicon cage. Accordingly we labelled the results corresponding to hexagon geometry with '*' 
and otherwise mentioned the discussions are for TM atoms that are bound in cage geometry. 

TM atoms can settle inside H-SiNW(N) without deforming the wire structure for core sites 
(s=3-5) and no structural difference is observed for different facets ([3]). On the other hand, edge 
effects become more important for subsurface sites. While no significant structural modification 
is observed for H-SiNW(61)+TM system, the Si-Si bond is broken upon TM impurity for H- 
SiNW(45)+TM(l) and H-SiNW(57)+TM(l) as shown in [3} The minimum and maximum dsi-xM 
ranges between 2.2-2.9 A(generally 2.4-2.8 A) and are summarized in[T]and[2[ for subsurface and 
core sites, respectively. 

The comparison of Ex indicates that the bond-breaking at site(l) for H-SiNW(45) and H- 
SiNW(57) results in an increase in Ej and makes it less favorable to other possible subsurface site 
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(s=2) where there is no deformation. On the other hand for the case of H-SiNW(61), site(l) which 
is inside a sharp corner (edge), is the lowest energy configuration. For core sites, Ex's slightly 
vary while moving from subsurface to the center (3— J-S) and the influence of surface becomes less 
significant. 

We calculate the binding energy(Eb) of TM atoms for different sites by using the expression: 

Eb =ET[H-SiNW(N)]+E[TM]-ET[H-SiNW+TM] (1) 

in terms of the total energy of optimized H-SiNW(N) and H-SiNW(N)+TM and the total energy of 
the string of TM atoms having the same lattice parameter as H-SiNW(N)+TM, all calculated in the 
same supercell. Interestingly, the variation of Eb with type of TM atoms follows Friedel modeP 
for both different sites and facets Q J ^^H^ I Accordingly, the lowest Eb is obtained for Cr and Mn and 
the highest Eb is obtained for Co and Ti depending on the number of filled d-states. When different 
facets are compared, Eb for N=61 is higher than N=57 and N=45 for subsurfaces sites, and there 
is no significant variation noticed for core sites. When doping sites are examined, interestingly 
the highest Eb is obtained at site(5), site(2), and site(l) for H-SiNW(45), H-SiNW(57), and H- 
SiNW(61), respectively. The results point out that the structural properties as well as the energetics 
depend on both the type of TM and also the cross section due to surface effects. 

Magnetic Properties 

The total energies(ET) for the considered systems are obtained both from spin-polarized (sp) and 
-unpolarized (su) states and the energy difference (AE = E^" — Ej) is used to determine the lowest 
energy configuration. According to this definition sp is the lowest energy configuration when 
AE > 0. As shown in [T] and |2} expect for Ti, the lowest energy configuration is obtained for sp 
state. Additionally we also compared the Ej's in ferromagnetic and anti-ferromagnetic states and 
confirmed that the ground state is ferromagnetic except Ti which yields paramagnetic ground state. 
The magnetic moments, jU vary depending on both the type of TM and also doping site as listed in 
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[Hand 11 

Using the energy difference between ferromagnetic and anti-ferromagnetic states, the Curie 
temperature (Tc) can be roughly estimated in the mean field approximation. Accordingly we con- 
sider H-SiNW(45)-i-TM as prototype, and calculated Tc as ~600K ,900K , and 400K, for Cr, Mn, 
and Fe respectively which indicate that ferromagnetic state is stable over room temperature. On the 
other hand Tc for V and Co is around ~50K which shows that they can only have ferromagnetic 
ground state at very low temperatures. Finally Ti (and also Ni) yields paramagnetic ground state 
even at OK. When the trend of Tc is compared with the unpaired J-electrons (d"") of TM dopants, 
a strong relation can be noticed. When d""=2 (Ti and Ni), Tc=0 (in other words paramagnetic), 
d""=3 (V, Co) Tc ~ 50, J""=4 (Fe) Tc ~ 400 and d""=5 (Cr and Mn), Tc > 600. Accordingly 
we can conclude that the stability of ferromagnetic state is mainly determined by the number of 
unpaired <i-states. 
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Figure 4: (Color online) The binding energy trend with respect to type of transition metal atoms 
(Ti, V, Cr, Mn, Fe, and Co) in H-SiNW(N) for N=45, 57, and 61 for the sites 1, 2(subsurface), and 
5(core). 



Electronic Structure 

Our calculations reveal that electronic structures are strongly affected by the type of TM atoms, 
doping site and also cross section of H-SiNW(N). The resulting band structures for subsurface 
(s=l,2) and core sites (s=4,5) are shown in Fi^5]j8] Depending on the configuration the ground 
state can be paramagnetic or ferromagnetic and electronic structure can be metal, semiconductor, 
semi-metal or half-metal.'^^^ 
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If we analyze the systems one by one, H-SiNW(N)+Ti configuration is paramagnetic metal for 
all considered systems. Interestingly, for V at site(l), the nanowires become half- metal, metal and 
semi-metal for N=45, 57, and 61, respectively and this indicates the effect of cross section on the 
electronic structure. Half-metallic state is also noticed for N=57 at site(2), and for all other doping 
sites (s=3-5), the system is ferromagnetic metal with varying magnetic moment. Cr and Mn behave 
very similarly. The nanowires are in general ferromagnetic metal for subsurface sites and they 
become half-metallic for core regions. Exceptionally H-SiNW(61)-i-Cr(4), H-SiNW(61)-i-Mn(4,5) 
are ferromagnetic semiconductor as TM is bound to hexagonal site instead of cage configuration 
(|3]). For Fe, all H-SiNWs are ferromagnetic semi-metal at site(l), but they become ferromagnetic 
semiconductor for all other sites including hexagonal doping site. And finally, for Co, all H-SiNWs 
are paramagnetic metal at site(l) and become half-metal or ferromagnetic semiconductor for core 
regions. 

Combining the results from the previous sections we can conclude that: 
(i)The energetics of binding is almost identical at core sites for different types H-SiNWs but the 
variation is noticed for subsurface sites due to surface effects (ii)The trend of Eb follows Friedel 
model^^ independent of H-SiNW type and binding site. (iii)The electronic ground state of H- 
SiNW(N)-i-TM is mainly determined by the type of (or the d-electron configuration of) TM dopants 
but can also differ for the same TM depending on the doping site and geometry. (iv)TM doping 
generally induces metallization (either paramagnetic or ferromagnetic) except for the case of Fe. 
H-SiNW(N)s still remain semiconducting after Fe doping with varying electronic band gap. The 
semi-metallic behavior of H-SiNW(N)-i-Fe(l) is an exceptional case and will be further discussed 
in the next section, (v) H-SiNW(N)-i-(Cr, Mn) systems are in general ferromagnetic metal for 
subsurface sites and they start to posses novel half-metallic behavior for core sites. (vi)However, 
H-SiNW(N)-i-V and H-SiNW(N)-i-Co also have half-metallic ground state for specific cases, they 
are ferromagnetic only at very low temperatures which make them practically paramagnetic metals. 
Accordingly, stable half-metallic ground state is unique to Cr and Mn where all the d-electrons are 
unpaired. vii)When Mn and Cr is settled in the center of silicon hexagonal plane, dsi-xM becomes 
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equal for all six nearest silicon atoms and it leads to different hybridization of p- (that belongs 
to H-SiNW) and d-orbitals (that belongs to TM) which makes the systems semiconductor instead 
of half-metal. However this configuration is energetically less favorable when compared to cage 
geometry but can be stable. 



N=45 N=57 N=61 N=45 N=57 N=61 N=45 N=57 N=61 




Figure 5: (Color online) The electronic band structures of H-SiNW(N)+TM at site(l) for N=45,57, 
and 61 and TM=Ti, V, Cr, Mn, Fe, and Co. Solid and dotted lines indicate spin up and down states, 
respectively. 



Hybrid functionals and GGA+U 

It is well-known that the usual approximations to DFT, such as LDA or GGA, lead to signifi- 
cant underestimation of the band gaps due to intrinsic failure in handling self-interactions.'^ GW 
corrections'^ appeared as a powerful tool to overcome this problem and were recently applied to 
silicon nano wires. f^^^ Unfortunately such GW calculations are computationally very costly and it 
is difficult to apply this technique on large systems. Moreover, GW corrections are performed on 
the frozen structure and do not allow for self-consistent structural relaxations. 

Alternatively, calculations using hybrid functionals,'^ which basically combine Hartree-Fock 
and DFT, usually provide improved electronic properties compared to DFT, although they are often 
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N=45 N=57 N=61 N=45 N=57 N=61 N=45 N=57 N=61 




Figure 6: (Color online) The electronic band structures of H-SiNW(N)+TM at site(2) for N=45,57, 
and 61 and TM=Ti, V, Cr, Mn, Fe, and Co. Solid and dotted lines indicate spin up and down states, 
respectively. 



N=45 N=57 N=61 N=45 N=57 N=61 N=45 N=57 N=61 




Figure 7: (Color online) The electronic band structures of H-SiNW(N)+TM at site(4) for N=45,57, 
and 61 and TM=Ti, V, Cr, Mn, Fe, and Co. Solid and dotted lines indicate spin up and down states, 
respectively. 
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N=45 N=57 N=61 N=45 N=57 N=61 N=45 N=57 N=61 




Figure 8: (Color online) The electronic band structures of H-SiNW(N)+TM at site(5) for N=45,57, 
and 61 and TM=Ti, V, Cr, Mn, Fe, and Co. Solid and dotted lines indicate spin up and down states, 
respectively. 

less accurate for structural properties. Hybrid calculations are less costly than GW corrections and 
also allow for structural relaxations. Hybrid functionals were used to calculate the electronic struc- 
ture of carbon nanotubesl^^^ and graphene nanoribbonSjI^and yield good agreement with exper- 
imental data. More recently, Rurali et al.'^ also applied hybrid functionals to silicon nanowires to 
accurately determine electronic band gap for varying diameter and orientation. 

In this section, we compared the results obtained in the previous sections with those obtained 
by using Bl-WC hybrid functional. We applied this technique only to specific cases, such as H- 
SiNW(45)-i-Fe at site(l) which is a ferromagnetic semimetal, H-SiNW(45)-i-Fe at site(5) which is 
a ferromagnetic semiconductor, and H-SiNW(45)-i-Cr [which is very similar to H-SiNW(45)-i-Mn] 
at site(5) which is a half-metal and compare the results obtained with GGA and GGA-i-U'^ calcu- 
lations as shown in|9l 

When Bl-WC is used, the band structure of H-SiNW(45)-i-Fe at site(l) is altered and becomes 
semiconductor instead of semi-metal. Similar band profile with a smaller band gap is also obtained 
when Hubbard-U (U=3 eV) term is considered for the strong on-site 3d electronDelectron interac- 
tions on Fe. Similar results are obtained for H-SiNW(45)-i-Fe at site(5). The system remains to be 
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Figure 9: (Color online) The electronic bands structures of H-SiNW(45)+Fe at site(l), H 
SiNW(45)+Fe at site(5), and H-SiNW(45)+Cr at site(5) obtained by Bl-WC and GGA+U cal 
culations. Solid and dotted lines indicate spin up and down states, respectively. 
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semiconducting but bang gap significantly increases when compared to GGA results. Accordingly 
we conclude that semi-metallic behavior of H-SiNW(N) is an artifact of DFT at GGA level and 
H-SiNW remains semiconducting with a modified band gap upon Fe doping for both subsurface 
and core sites. 

The novel half-metallic ground state of H-SiNW(45)-i-Cr at site(5) is preserved by B 1-WC and 
also when U term is introduced. Interestingly, the dispersive metallic spin state (f) is not affected 
by either approaches but the the band gap for the insulating spin state (J,) is modified. In other 
words, half-metallic ground states obtained at the DFT-GGA level™ are even more stable than 
predicted. 

Finally, the results indicate that Bl-WC is not only capable of correcting the under estimated 

electronic band gaps but also can handle strong on-site d electron-electron interactions without 

introducing Hubbard-U term, being a good alternative to LDA(GGA)-i-U for strongly correlated 
systems.lMI 

Conclusions 

We analyzed the structural, electronic, and magnetic properties of hydrogen saturated silicon 
nanowires with interstitial transition metal doping. We found that the electronic and magnetic 
ground state is mainly determined by the type of (or d-electron configuration of) transition metal 
atom, it can also be affected by doping site and cross section (surface effects) of silicon nanowire. 
Upon transition metal doping there is a tendency for metallization except for the case of Fe and 
even novel half-metallic configuration can be obtained for specific cases. The stability of ferromag- 
netic ground state and energetics of binding are inversely proportional and depends on the number 
of unpaired 3J-electrons following the Friedel model. The obtained results are also compared with 
those calculated using new type hybrid functional (Bl-WC) in order clarify the limitations of den- 
sity functional theory at GGA level and also to find out the capabilities of Bl-WC. We believe that 
our detailed analysis will guide both experimental and theoretical studies related with doping of 
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silicon nanowires and other systems and moreover present results hold the promise for the use of 
silicon nanowires in various spintronic applications upon transition metal doping. 
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